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We study the large-time behavior of a class of periodically driven macroscopic systems. We find, 
for a certain range of the parameters of either the system or the driving fields, the time-averaged 
asymptotic behavior effectively is that of certain other equilibrium systems. We then illustrate with 
a few examples how the conventional knowledge of the equilibrium systems can be made use in 
choosing the driving fields to engineer new phases and to induce new phase transitions. 
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^ ■ I. INTRODUCTION 

o 

With an enduring interest studies have been pursued to understand systems with many degrees of freedom at, 
I i near, and far away from equilibrium. Systems at equilibrium are studied within the well-established Boltzmann- 
Gibbs framework, while those near and far from equilibrium lack such a framework to study within and are usually 
described by stochastic dynamical equations. The stochastic equations that describe the dynamics near equilibrium 
need to ensure that the system relaxes to equilibrium and hence are less difficult to construct than when the system is 
d . far away from equilibrium. However this condition, that the asymptotic solution is the equilibrium distribution, does 
not provide a unique stochastic equation. In fact, for a given equilibrium system infinitely many such equations near 
I . equilibrium can be provided that would lead to the same static properties [1] . This naturally motivates the study of 
various stochastic models. 

We know that a closed system, characterized by the Hamiltonian H, when exposed to the environment for a long 
time, will equilibriate and be described by the Boltzmann distribution exp(—f3H). If this system is also subjected to 
an external periodic force then what is its large-time behavior? Due to the absence of an established framework to 
I ' find this asymptotic behavior it is essential to study various periodically driven stochastic systems. It is also essential 
, as some of these systems show interesting behavior and some of them furnish useful applications. Stochastic resonance 
C ' [2] and magnetic hysteresis [3] are some of the well known phenomena that are exhibited by certain periodically driven 
stochastic systems. 

A large number of stochastic processes have been studied for many decades now [4]. Studies in critical dynamics 
have been actively pursued [5] after renormalization group techniques were successfully applied to equilibrium systems 
[6]. Many driven diffusive models were studied [7] and the critical behavior of some of them were analyzed too [8,9]. 
Periodically driven stochastic particle systems were extensively studied in the context of stochastic resonance [2], 
while similar studies on fields were relatively few. 

The issues that we address in this paper are related to a class of periodically driven macroscopic systems that, in 
the absence of driving, would relax to equilibrium. The relevant degrees are described by a field and the equilibrium 
, properties are characterized by an energy functional of the field. These systems may or may not have a critical 
point when in equilibrium but could exhibit a critical behavior when the driving fields are switched on. Some of the 
questions that we ask: How do these macroscopic systems respond to periodically driven fields? How do we describe 
the phases of the driven system? What is the dependence of the critical points, if any, on the parameters of the 
driving fields? Can driving change the nature of the phase transitions? Can it lead to new phases? 

The layout of the paper is as follows. In the next section, we define a class of stochastic models to describe the 
relevant fields of periodically driven systems. In Sec. Ill, we develop a perturbative scheme to solve the Fokker-Planck 
equation that describes the systems of our interest and explicitly find their asymptotic behavior to first order. In 
Sec. IV, we illustrate with examples some of the effects of driving on phases and phase transitions. 



II. STOCHASTIC MODELS FOR PERIODICALLY DRIVEN SYSTEMS 

In this section, we will define a class of periodically driven stochastic systems. We will assume the slowly relaxing 
modes of the system in the absence of the periodic forces to be the relevant degrees of freedom and then model the 
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dynamics of these variables by modifying the stochastic equations of these modes in the presence of the periodic 
forces. 

Let us first recollect how the dynamics near equilibrium is described in the absence of driving. In this case the 
variables of interest, those that relax very slowly, are the order parameter and the conserved quantities, if any. This 
is a small set of variables which are some functions of the original degrees of freedom. In principle, the dynamics of 
these relevant variables can be obtained from the equations of motion of the original degrees, by integrating out the 
unwanted variables. In practice, these dynamical equations for the relevant degrees are not thus established as they 
are not usually derivable due to the large number of unwanted variables involved. Hence these equations are based 
on certain guiding principles [10]: (i) Forgoing a large number of unwanted variables renders a stochastic evolution to 
the relevant ones, (ii) These stochastic dynamical equations should be such that the large-time distribution gives the 
right static thermodynamic properties. This large-time distribution is the distribution either obtained by integrating 
out the irrelevant variables from the canonical equilibrium distribution or found on phenomenological grounds, for 
e.g., Landau- Ginzburg (LG) theory near critical point. Finally, the stochastic equations thus established could explain 
the dynamic properties of the relevant variables near equilibrium. 

There exists many different families of stochastic equations that give the same static properties for the relevant 
degrees. One of these is the following Langevin equation, the time-dependent Landau-Ginzburg (TDLG) theory that 
describes the dynamics near equilibrium, 



r|,(M)- SnM 



dt ' Sip (a 



+ V(x,t), (1) 



where tp(x, t) is a d + 1 dimensional stochastic field that relaxes to the order parameter ip(x) of the system and 
H(<p) is the LG free-energy functional; rj(x,t) is a Gaussian random field with (77(2;, i)),, = and (t](x,t)r](x' ,t'))r) = 
2r<5(a; — x')S(t — t'). If p v (x,t)) is the solution of the above Langevin equation then the probability distribution of 
this solution, P((p,t) = Y\ x (S(<p(x) — ip v (x,t))) ni evolves according to a Fokker-Planck (FP) equation, d t P = CP, 
that is obtained from Eq.(l) and the Gaussian distribution of r](x,t). The formal solution of this FP equation is 
P((p,t) = cxp(Ct)P(ip,t = 0), where P(<p>,t = 0) is the initial distribution. The asymptotic distribution, Poo(ip,t) = 
lim^oo P(ip, t), is the right eigenfunctions of L operator with zero eigenvalue. This asymptotic probability distribution 
is the LG measure, P eq (<p) ~ exp(— H(ip)) , provided (i) the eigenvalues of C are negative semi-definite, (ii) there exists 
a unique normalizablc eigenfunction corresponding to zero eigenvalue in the connected space of functions to which 
the initial distribution belongs to, and (Hi) the initial distribution is not orthogonal to this eigenfunction. 

The presence of driving fields on the system will alter both its kinematics and dynamics. Assumptions about 
kinematics are: (i) The set of variables that were relevant in the absence of driving for the dynamics near equilibrium 
continue to remain relevant upon introducing driving, (ii) Some minimal set of additional variables might become 
relevant too. We will consider a velocity field n(x) along with the order parameter field ip(x) to be relevant. The 
dynamical variables of the system are then specified by ip(x,t) and its time derivative dt<p(x,t). Assumptions about 
the dynamics are: (i) The Langevin equation, that describes dynamics near equilibrium in the absence of driving, gets 
modified by adding terms related to the periodically driven fields and by adding a second-order time-derivate term 
dfip(x,t) whose effect is significant when the driving frequency is high, (ii) The periodic driving does not change the 
properties of the noise field for any frequency of the driving fields. 

The modified TDLG theory that describes the system when subjected to driving is then given by the following 
stochastic equation 

m^(x, t) = -Y ^<p(x, t) - + F{<p{x, t), t) + V (x, t), (2) 

where !F(ip(x), t) is the driving field that is periodic in time with period T = 2n/Cl. The general form of this field is 



F(v(x),t)=Y J 



71=1 



F n (ip(x)) cos(rafit) + G n (ip(x)) sin(n£li) 



(3) 



where F n and G n are arbitrary local functionals of <p(x). This term could be thought of as a result of making the 
coupling constants in H(<p) time dependent and periodic. We will take m — 1 so that the expressions look simple and 
to regain the m dependence back replace t — > t/^/m, T — > Y / y/m and Q \/m£l. The distribution of the noise field 
is Gaussian as specified earlier. 

The above dynamical equation for a 0+1 dimensional stochastic field describes periodically driven Brownian particle. 
These Brownian particles exhibit a variety of interesting asymptotic behavior depending on the driving forces [11]: 
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The particles when driven could congregate around more than one point even though when in equilibrium they would 
have congregated around a single point. Particles with different masses respond to driving differently and thus when 
they are mixed and driven would cluster around different points. 

The phase space probability distribution is defined as P((p,Tr,t) = (Y[ x S((f(x) — <p^(x, t))<5(7r(x) — dt(p v (x, t))) , 
where ip^x^t) is the solution of Eq.(2) at time t for a particular history of {rj(x,t)} over a time t and <■■■>,, is 
the average over the noise distribution. This distribution satisfies the normalization condition J T>{ip, ir]P((p, n,t) = 1. 
The time evolution of P(ip, n, t) is described by the FP equation 

d -P(ip,TT,t) = C(t)P(<p,ir,t), (4) 



where 



<(t) = I 

J X 



dt 



5 6 / SH r . A „ 5 2 



= C 



FP 



Li™- (5) 

The FP operator Cfp is defined as the time- independent part of the C(t) operator, both of which are assumed to be 
suitably regularized. To keep the notation compact, <p — <p(x) and it = ir(x) are used in the above equation. The 
solution of this equation is not known in general and our aim is to find the asymptotic solution for some range of 
parameters of the system and driving fields. 

The asymptotic distribution of the above FP equation is a periodic function with period T, if the real part of 
the eigenvalues of dt — C(t) are positive definite. In brief, the argument for the periodicity goes as follows. The 
solution of the FP equation is the right eigenfunction of the operator dt — C(t) corresponding to the zero eigenvalue. If 
C{t) = C(t + T) then d t — C(t) commutes with the discrete time-translation operator exp(T<9 t ). The solution can then 
be expanded in terms of the common right eigcnfunctions of these two operators. Let the eigenfunctions of dt — C(t) 
and exp (T<9 t ) with eigenvalues and exp (/xT), respectively, be the Floquet-type functions exp (/rf)p^(i), where 
is a periodic function with period T. Substituting these eigenfunctions in FP equation gives [d t — C(t)]p^(t) = —fip^t). 
Hence, if the real part of the eigenvalues of dt — C(t) are positive definite then in the limit t — > oo the only eigenfunction 
that survives is po(t), thus making the asymptotic distribution periodic. 

When the coupling constants {g} of the driving field !F(ip,t) are small compared to the coupling constants {g} of 
Ti.(ip) then the FP equation can be solved perturbative in {g} provided we know the right and left eigenfunctions of 
the FP operator Cfp- These eigenfunctions are generically not known though the eigenfunctions of the FP operator 
which includes only the free part of Tt((p) are obtainable. Hence, the eigenfunctions of Cfp can be determined 
perturbatively in {g} and in turn the solution of the FP equation in the double series expansion in {g} and {g}. 

We will now further to find the asymptotic behavior of the FP equation when {g} are not necessarily small compared 
to{ 3 }. 



III. ASYMPTOTIC DISTRIBUTION OF THE FP EQUATION 

In this section, we will first transform the FP equation to enable us to have a non-perturbative solution in {g} 
but perturbative in a parameter that involves both {g} and £l 2 + T 2 . We then explicitly evaluate the asymptotic 
distribution to first order in this parameter and more specifically the time-averaged asymptotic correlation functions. 
We find that these correlation functions can be expressed as equilibrium correlation functions with an effective LG 
energy functional. 



A. Formal solution 



We will make a change of variables and transform the FP equation such that the time-dependent part of the 
transformed FP operator is small though the time-dependent part of the original FP equation is not. The variables 
{ip, it, t} are changed to {$, n, r} by the following nonlinear transformation 

<p(x) = $(x)+£(x;&,t), 
ir(x) = U(x) + d T £(x; $, r), 

t = r, (6) 
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where the explicit form of £(x; $,t) will be specified later. Under this transformation the probability distribution is 
made to behave like a scalar: P(tp, n, t) — > IT, r) = P(<£, 7r, t). The functional derivatives will then transform as 



5ir(x) 5Ii{x) ' 
5 

5ip(x) 



M{y,z) 



5U(z) J ' 



and the time derivative as 



where 



a_ _ d_ 

di~dr 



M{x,y) 



S(p(x) 



5(x — z). 



(7) 



(8) 



(9) 



Substituting the above new variables in Eq.(4) we obtain the following transformed FP equation 



dP 
~dt 



-L 

+ [ n(x) 

Jx.y 



-4 n - Jn (- rn + ^ + A + A ^< r )) + r ^ 



8<f> SU 

s ( x -y)~ D{x,y) 



S 



d*n + rd T £-f($,T) 



5TL(x) 



S$(y) 
P. 



D(x,y) I 

J z 



5IP 



d T M(y,z) 



8U(z) 



P 



(10) 



where fn(f( x )) = —5H/5(p(x) and A.F(£, r) = JF($ + £, t) — T(<&, r). We now choose $, r) to be the solution of 
the following equation 



whose explicit form is 



-1 



n 2 tl 2 + r 2 



^ + ra T e-^(*,r) = o, 



(*(*)) + ^G„(*(a;)) ) cos(nOT) + [G n ($(x)) - —F n ($(x)) ) sin(nf>T) 



(11) 



• (12) 



The reason for this specific choice for £ is that not only does it make the last term on the right hand side of Eq.(10) 
zero but also makes the (time)r-dependent part of the modified FP operator to be of 0(£). Thus if £ is small then 
this modified equation will yield a perturbative solution. Upon substituting £ from Eq.(12) in Eq.(10) we finally get 
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&r -p($,n,r) = [£ + A£]P(<f>,n,T), 
where C is the following static FP operator 

£ - LH n - m (-™ +/^>+^>) + *m 

and A£ is the periodic (time)r-dependent operator 

AC = J x -jri ( fn ^ + -/«(* + + - rj) 



(13) 



(14) 



+ ^ n(ar) ^8(x-y)-D(x,y)^-J^+D(x,y) J^d T M(y,z) 



SIL(z) 



(15) 



where the bar over the terms indicate an average over a time period. 
The perturbative asymptotic solution can be formally written as 
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P 00 ($,n,r) = Q($,n) + ^-^A£P 00 ($,n,r), (16) 

where Q(Q, II) is the right eigenfunction of L with zero eigenvalue. If L has a unique cigenfunction corresponding to 
zero eigenvalue, and the_real part of the non-zero eigenvalues do not vanish, and P M is periodic in time, then it follows 
from Eq.(13) that ACP^ has no overlap with the eigenfunction of d T — C corresponding to zero eigenvalue. Thus, 
though d T — C is not invertible, its inverse action on the space orthogonal to its eigenfunction with zero eigenvalue is 
well defined. 



B. Effective theory 

In this subsection, we will obtain the effective energy functional to 0(£) by averaging the asymptotic correlation 
functions over a time-period. The observables will be related to the time-averaged correlation functions if the time of 
observation is comparable to the time period of the driving fields. In other words, we are assuming that the driving 
fields oscillate rapidly compared to the time scale of the measurement. 

The equal-time correlation functions of the stochastic field transform under the change of variables as follows 

(ip(x 1 ,t)ip(x 2 ,t) ■ ■ -) n = J D[tt, ip]<p(x 1 )ip(x 2 ) ■ ■ ■ P(ip, it, t) 

= J P[n,$]J[$]|$( Sl ) + ^ Xi;(i >,i)||$( S2 ) +e ( X2;(i >,i)|...p($,n,i), (17) 

where the integration measure T>[Tr,(p] — Y[ x [d 7T (x)d(p(x)] and the Jacobian of the transformation J[<f>] = Ylxi-^ + 
d^(x, <&(x), t)/d$(x)]. The first-order asymptotic distribution, obtained from Eq.(16), is Poo = + (d T — 

C)~ x A£Q(°\ where and are the solutions of the equation CQ — to O(l) and 0(£), respectively. Substi- 
tuting Pqo in Eq.(17) and averaging over a period gives 

J %,*(nMi 2 )-ftSM = / ^[n,$]$(x 1 )$(x 2 )--.g( 1 )($,n) + o(e 2 ), (is) 



since the time-average of both £ and [(d T - C^ACQ^} are zero. From Eq.(14), in which fn{$ + and AJF(£,t) 
are expanded to 0(£), we get 

g( 1 )(d>,n) = ^ TT cxp(^-i |n 2 (,)-H($)-AH($)), (19) 

where is the normalization constant and AH(Q) satisfies the condition 



AW($) = -^;$,t)-A-f($(a : ),t). (20) 



S$(x) v ' sv ' ' 'd$(x) 
Substituting £ and T from Eqs. (12) and (3) in the above equation we get AH = AHi + AH 2 , where 



= \ E n 2 Q 2 +T 2 I (*(*)) + G l 

n=l Jx \ , 

^ (*) - E {n W + T*)(2r*l) L W ^ ' G " ] (21) 

W[F n ($(x)),G n (<f>(x))] is defined by the indefinite integral W[F(a),G(a)] = J da[G(a)F'(a) - F(a)G'(a)}. After 
integrating out II in Eq.(18), we finally get 

lim {ip{x u t) V {x 2 ,t)---)r, = \\ T>[ip{x)]<p(x x )<p{x2) ■ ■ ■ e- H °f'M, (22) 

t^oo Z J 

where Z = fT>ipexp(— H e ff) and H e //(^) = H((p)+AH(<p). The effective energy functional for a TV-component order 
parameter field ip = {<p a } is a straight-forward generalization of Eqs. (20) and (21). In this case AH, which satisfies 
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Eq.(20) for each component, exists only if the components of the driving field F% and G° satisfy the conditions, 
dFZ/dtpt, = dF b Jdy a and dG a Jd Vb = dG b Jdip a . 

Note that we have neglected 0(£ 2 ) terms in the expansion of the C and A£ operators because £ is assumed to be 
small. The condition for the smallness of £ and the criteria for the validity of 0(£) approximation arc both obtained 
by comparing the neglected terms with those that are retained in these operators. For instance, if T — 2gip 2 (x)cos(flt) 
and the coefficient of p 4 term in H(p) is A then £ is small if A 3> g 2 /(£l 2 + T 2 ). A self-consistent criteria can also be 
obtained by comparing the expectation values of the neglected and the retained terms using Ti e f f . In some cases the 
higher powers of £ contribute only irrelevant terms to Ti. e f / and hence can be neglected while determining universal 
properties, independent of £ being small. For instance, if T = g(p 4 (x)cos(Slt) and H is a ip 4 theory in four dimensions, 
then the 0(£) term in 7i e // itself is an irrelevant p 8 field and the higher-order terms contain fields that are more 
irrelevant. 

The effective energy functional Tt e ff can be interpreted as follows. The asymptotic distribution contains both 
statistical fluctuations $ and the dynamical fluctuations t) that are periodic in time. In Eq.(20), if we substitute 
T in terms of £ by using Eq.(ll) and then integrate by parts, we get 

Thus, by averaging over the time period we have eliminated the dynamical fluctuations and provided an effective 
description to the system in terms of the statistical fluctuations that are governed by a modified energy functional 
H eff =H + AH. 



IV. ILLUSTRATIVE EXAMPLES 

We have seen that the large-time behavior of periodically driven stochastic systems, averaged over a period, can 
be described by an effective LG functional of equilibrium systems. We can now make use of the knowledge about 
these equilibrium systems to induce new phases and phase transitions in various systems by subjecting them to 
time-dependent periodic fields. We illustrate some of the effects due to these driving fields with a few examples. 



1. Varying the critical point 

The driving fields transform the energy functional 7i into a new energy functional H. e f f . This amounts to changing 
the coupling constants, {g} — ► {g e }, which in turn induces a change in the critical point. We will illustrate this with 
an example. The LG functional that describes the Ising model near the critical temperature is 

H[p]= f(d<p(x)) 2 +A(T-T c )p 2 (x) + \p 4 (x). (24) 



The mean-field theory suggests a second order transition at T c from a Z 2 -symmetry-broken phase to the symmetry- 
unbroken phase and a dependence of magnetization (order parameter) on temperature below T c as (i^) 2 = A(T C — 
T)/2A. If we now drive this system, say, by oscillating ip 2 term, which is same as adding a force term 

F(<p(x)) =2aip(x)cos(Q,t), (25) 

then this system will get described by the following effective energy functional that is obtained using Eq.(21) 

We// M = I (M*)) 2 + aeV 2 (x) + \<p 4 (x) , (26) 

J X 

where a e = A(T — T c ) + a 2 /(0 2 + T 2 ). Since this effective energy functional differs from the original functional only 
in the coefficient of the ip 2 term one can read off the critical temperature 8 C and the behavior of the time-averaged 
magnetization (ip) below 9 C of the driven system. We get 

Or. = Tc. — 



A(Sl 2 +T 2 ) : 



(pf = ^(e c -T). (27) 

Hence, the driving field acting on this system tends to destroy the symmetric phase as it reduces both the critical 
temperature and the magnetization at a given temperature. 
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2. Changing the nature of transition 

The nature of phase transition can also be changed by applying driving fields. Consider the LG functional as given 
in the previous example, Eq.(24), but with a different driving force 

T(<p{x)) = 2(aip(x) - bip 3 {x)) cos(fii), (28) 

and also assume the dimension of space d > 3 where ip 6 term is irrelevant. The effective energy functional will then 
become 

W e //M = / {9<p(x)) 2 + a eV 2 {x) + \ e tp 4 (x) + b e <p 6 {x), (29) 

J X 

with the coupling constants a e = A(T - T c ) + d 2 /(Q 2 + T 2 ), A e = A - 2ab/(fl 2 + T 2 ), and b e = b 2 /(Q 2 + T 2 ). 

Now the mean-field prediction is as follows. For A e > two different phases exist; a spontaneously broken phase 
when a e is negative and an unbroken phase when a e is positive. These phases are separated by a line of second order 
transition which ends in a tricritical point when both A e and a e become zero. For A e < there are three possibilities 
depending on the value of a e : (i) If a e < \ 2 /4b e then it is in a symmetry broken phase, (ii) If A 2 /46 e < a e < A 2 /36 e 
then it is a metastable phase and will have a broken symmetry if one enters this region crossing the line A 2 = Aa e b e 
and will have an unbroken symmetry if one enters by crossing the line A 2 = ?>a e b e . (Hi) If a e > A 2 /36 e then it is in 
the symmetric phase. The line A 2 = Aa e b e is a line of first-order transition which ends in the tricritical point. 

Suppose the system initially is in the symmetric phase with both a e and A e positive. Now switch on the periodic 
force such that the product ab is positive and then gradually reduce the driving frequency f2. This will reduce both a e 
and A e and takes the system across the line of first-order transition into the symmetry broken phase. More generally, 
by tuning at most two of the three parameters {a, b, ft} and the temperature T we can scan the entire a e — A e plane. 
Hence, by applying the driving fields we can engineer the behavior of the system. 



3. Inducing new fixed points 

Driving fields can enlarge the coupling constant space of the system by introducing either relevant or irrelevant 
fields. Though the irrelevant fields usually do not change the large-distance properties of the system but in some 
cases, like in the previous example, the coupling constant of the relevant fields can get shifted to a region where 
irrelevant fields become important. Let us now examine an example where the driving fields introduce a relevant field 
and drastically alter the large-distance properties. This is because the system flows under scaling to a new stable 
fixed point in the enlarged coupling constant space. The system under study is described by the O(N) model 

H = [ [d<p{x)f + t v 2 (x) + u(^ 2 {x))\ (30) 

J X 

where ip = {ipi, • • • , lpn} is a iV-component vector field. This model has two fixed points: one is at u = (Gaussian 
fixed point) and the other is at finite u (Heisenberg fixed point). For dimension d < 4 Heisenberg fixed point is stable 
while Gaussian is not. Now drive the system by the fields 

T a ( V (x)) =2^ 2 (a;)cos(^), (31) 
a = 1, • • • , N. The resultant model is the 0(N) model with a cubic symmetry breaking term [12]: 

H eff =H + v f^vHx), (32) 

where v = (f2 2 + r 2 ) -1 . In the u — v plane this model has two more fixed points, Ising fixed point at (u = 0,v ^ 0) 
and the cubic fixed point at (u ^ 0,!) 5/ 0), apart from the Gaussian and Heisenberg fixed points at (u = 0,v = 0) 
and (u ^ 0,v = 0), respectively. The stability of the fixed points is as follows: The Gaussian point is unstable in 
the ^-direction too. Ising point is stable in the v-direction but unstable in the w-direction. Heisenberg point is stable 
in the w-direction if N < 4 and unstable if TV > 4. Cubic point is stable in both directions if Heisenberg point is 
unstable and vice-versa. Thus, when TV > 4 these driving fields change the large-distance properties of the system 
from Heisenberg to cubic. Periodically driven stochastic models with O(N) symmetry, without the md 2 <p(x,t) in 
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Eq.(2), were first studied in the context of magnetic hysteresis [13,14]. The phase transitions that get induced by the 
driving fields in these models were recently investigated too [15]. 

In summary, we have derived the effective theory for the correlation functions of a class of periodically driven 
macroscopic systems. We have shown with a few examples that this effective theory can be made use to select the 
driving fields that can steer the system through a plethora of phases. 
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